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We present a search for the lepton flavor violating decay B + — > K + t^ // using 383 million BB 
events collected by the BABAR experiment. The branching fraction for this decay can be substantially 
enhanced in new physics models. The kinematics of the tau from the signal B decay are inferred from 
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the K + , n, and other B in the event, which is fully reconstructed in one of a variety of hadronic decay 
modes, allowing the signal B candidate to be fully reconstructed. We observe no excess of events 
over the expected background and set a limit of B(B + — » K + th) < 7.7 x 10~ 5 at 90% confidence 
level, where the branching fraction is for the sum of the K + t~ fi + and K + t + /j,~ final states. We 
use this result to improve a model-independent bound on the energy scale of flavor-changing new 
physics. 

PACS numbers: 13.25.Hw, 14.40.Nd 



Lepton flavor violation (LFV) thus far has only been 
observed in the neutrino sector P, 0, Manifesta- 
tions of LFV in S meson decays that have final states 
with charged leptons (e.g., S — ► KW) are allowed in 
standard model interactions, if massive neutrinos are in- 
cluded but such processes occur only at the one- 
loop level and are extremely suppressed by powers of 
ml/Myf- Branching fractions for lepton flavor violating 
B decays can be substantially enhanced in many exten- 
sions of the standard model [1, [E 0, S| • The semileptonic 
decay B — > Kt/i is likely to have higher sensitivity to 
new physics, when compared to leptonic decays such as 
B° — > t£' , since the latter is both helicity and Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix (CKM) sup- 
pressed by a factor of |Vtd/T4&| 2 . Some new physics mod- 
els require flavor changing neutral currents (FCNC's) to 
occur at the one-loop level, as in the standard model. In 
other extensions, such as models with a Z or additional 
Higgs doublets, FCNC's occur naturally at the tree level, 
unless they are eliminated by imposing an ad hoc discrete 
symmetry. 

A limit on the process B — > Krfi, which involves the 
second and third generations of both quarks and leptons, 
would provide a unique and powerful constraint on model 
parameters of grand unified theories. Cheng, Sher, and 
Yuan [Bl, 0] propose that in models with an extended 
Higgs sector the most natural value of the FCNC Yukawa 
couplings connecting generations i and j are proportional 
to ^ymiirij /m T , which implies that FCNC's in these the- 
ories should be largest in processes involving the second 
and third generations. An observation of B — ► Krfi 
would be an unambiguous sign of physics beyond the 
standard model. In this paper, we present the results of 
a search for B — > Kt[i. 

We use a data sample of 383 million SS pairs pro- 
duced by the PEP-II asymmetric-energy e + e~ collider, 
running at the T(AS) resonance, collected by the BABAR 
experiment [lj| at SLAC. Charged particles are identi- 
fied using a Cerenkov radiation detector and dE/dx mea- 
surements in the tracking system. Instrumentation em- 
bedded within the iron of the flux return for the 1.5 T 
solenoid aids in the identification of muons. An electro- 
magnetic CsI(Tl) crystal calorimeter (EMC) is used to 
reconstruct photons and identify electrons. 

The analysis strategy is to reconstruct the Y(4S) — > 
B+B- in the search for B+ -> K+TfJ, 11]. One of the B 
mesons (Stag) is fully reconstructed in one of a large num- 



ber of hadronic final states, B- -> The X~ 

represents a system of charged and neutral hadrons with 
total charge —1 composed of ni7r ± , n^K^, n^Kg, and 
r&4 7T°, with ri\ + n% < 5, rz.3 < 2, and n.4 < 2. The D*° is 
reconstructed in the D°7r° and S°7 channels, the D° in 
the K~ir + , K~ir + ir°, K~ir + ir~ir + , and K s 7r + ir~ chan- 
nels, and Kg in the 7r + 7r~ channel. We search for the 
decay B + — > K + Tfi using the remaining tracks in the 
event. The momentum vector of the signal B candidate, 
Psig, must be equal in magnitude and opposite in direc- 
tion to that of S tag in the center-of-mass (CM) frame. 
The t candidate kinematic variables, E T andp r , are fully 
constrained by p s j g , the measured momenta of the K + 
and [i tracks, and the constraint E T = -©beam — Ex — S„, 
where -Ebeam is the CM beam energy. The reconstructed 
r invariant mass m T = y/ E% — peaks sharply at the 
true r mass for the signal. 

Events are required to contain a Stag candidate with 

to es = \J £wm Ptag > 5 - 27 GeV/c 2 and Stag consis- 
tent with -Ebeam within three standard deviations, where 
Ptag and Stag are the momentum and energy of the S tag 
candidate in the CM frame. A S tag meson is fully recon- 
structed in about 0.2% of Monte Carlo events where one 
B ± decays to K ± Tfi. Even though the r daughters are 
not needed to reconstruct m T , we require the r in the 
signal B candidate to be consistent with a "one-prong" 
(i.e., one-charged-track) r decay to reject combinatoric 
background. Therefore we require exactly three charged 
tracks in the event not associated with S tag and with 
net charge opposite to that of S tag . Among these three 
tracks, we require: a kaon candidate with charge opposite 
St a g, a muon candidate, and a third track (the r daugh- 
ter) with charge opposite the muon candidate. The event 
is rejected if any of the three tracks is consistent with a 
proton hypothesis or if either of the two non-kaon tracks 
is consistent with a kaon hypothesis. Signal candidates 
are divided into three categories based on the properties 
of the r-daughter track: electron, muon, and pion. 

The kaon, muon, and electron particle identifica- 
tion criteria used in this analysis have momentum- 
dependent efficiencies and misidentification probabilities 
(fake rates). The kaon candidate must pass loose selec- 
tion criteria, based on the measured Cerenkov angle and 
dE/dx in the tracking system. Muon candidates, cither 
from the S decay or from the r decay, must pass mini- 
mum selection criteria that are 85% efficient for muons 
above 1.5 GeV/c and less than 10% efficient for pions and 
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kaons. Tau daughter electrons must pass minimum elec- 
tron selection criteria that are 95% efficient for electrons. 
More stringent electron and muon identification criteria 
for the tau daughter track are incorporated through a 
likelihood ratio described below. Tau daughters that do 
not pass either the electron or muon criteria fall into the 
pion category. 

For muon channel signal B candidates, there are two 
muons in the final state: one from the B decay (primary 
muon) and one from a leptonic r decay. Of the two pos- 
sible track assignments for the primary muon, we use the 
one that gives m T closest to the known r mass. The bias 
in the background m T distribution for the muon channel 
from using this procedure is found to be negligible. 

Semileptonic B decays can produce final states that ap- 
pear identical to the signal. For example B + — > D° 
followed by D° — > K + ir~ produces a K + h + -k~v final 
state. If the D° decays semileptonically, the final state 
is K + [i + £~ vv . These backgrounds are easily removed by 
requiring that the invariant mass, m(KTr), of the kaon 
candidate and the oppositely-charged, signal-track candi- 
date, when this track is assumed to be a pion, be greater 
than 1.95 GeV/c 2 . This requirement is greater than 50% 
efficient for the signal and removes about 99% of the 
background from BB events. 

The B+ -► and B+ -► D^fi+u^, with D° -> 

K + tt~ , backgrounds also form the D/iv control sample, 
which we use to normalize the signal branching fraction. 
Events for the D/iv control sample are required to have 
m(Kn) in the range [1.845, 1.885] GeV/c 2 (within about 
three standard deviations of the D° mass) . The neutrino 
momentum is calculated from p t ag and the three tracks 
in the D/jiv final state. We then compute 



AEu^v = Ek + +E-x+ Vv — Ebc&m — Vv — E n 



(1) 



We use AE DfJ/l , rather than m miss = v^miss Put similar 
to our m T reconstruction, because the expected D°fiv 
missing mass is zero. The AEd^v distribution for D°fj,v 
decays is centered at zero, while for D*°fiv events, it is 
shifted by —150 MeV and slightly asymmetric, due to 
the missing neutral particle from the D*° — > (7r°,7)_D° 
decay. We determine the yield of D° fxv and D*° fxv events 
simultaneously in an unbinned maximum likelihood fit of 
AE Dfn , (Figure [I]). 

For the K + t/j, signal, a non-negligible source of BB 
background in the sample remaining after the m{KTz) re- 
quirement comes from B + — ► (cc)K + with (cc) — ► 
where (cc) is a charmonium resonance. This background 
mainly enters in the muon channel, but is also present in 
the pion channel. For the muon and pion channels this 
background is removed by requiring the invariant mass 
of the two non-kaon tracks, when both are assumed to be 
muons, to be outside of the ranges [3.03, 3.14] GeV/c 2 
and [3.60,3.75] GeV/c 2 , which are centered on the J/tp 
and i()(2S) resonances masses respectively. 




-0 0.2 

AE D ^(GeV) 

FIG. 1: Results of the D[w control sample AEd^v fit. The 
points with error bars are the data, the solid curve is the 
projection of the fit, the dashed (dot-dashed) curve is the 
D°^v (D* fj,v) component, and the dotted curve represents 
events from other sources. 



At this stage in the selection, the background is dom- 
inated by continuum events (e + e~ — > qq with q = 
u,d,s,c). We suppress this background using a likeli- 
hood ratio 



Lr 



(2) 



where a? is a vector of four discriminating variables (see 
below) and P s (xi) (Pb{xi)) is the probability density 
function (PDF) , which describes the signal (background) 
for variable X{. The PDFs are separate for each of the 
three signal categories. The four discriminating variables 
are: | cos #thr | the magnitude of the cosine of the CM 
angle between the Stag thrust axis and the thrust axis 
of the rest of the event; ^2 E ca \ the total neutral EMC 
energy that is not associated with Btag5 the quality of 
the primary lepton identification; and the quality of the 
secondary lepton identification (for electron and muon 
channel candidates) . We require a minimum E c& \ energy 
of 50 MeV (100 MeV) for clusters in the barrel (forward 
endcap) to be included in ^E ca \. 

The | cos (9 t hr | distribution is flat for signal and peaks 
near one for continuum. The ^ E ca \ distribution peaks 
at zero for signal, while the background distribution is 
broad, peaking at around 1.5 GeV. The lepton quality 
is divided into four hierarchical, mutually-exclusive cate- 
gories with fake rates decreasing with increasing quality 
rank. For the highest-quality muon candidate rank, the 
fake rates from pions and kaons are less than 2%. The 
highest-quality electron candidate rank has a fake rate of 
less than 0.1% for pions, as high as 3% for low-momentum 
kaons, and below 0.4% for kaons above 0.8 GeV/c. We 
fit signal and background Monte Carlo histograms of 
| cos#thr| and ^cal to define P s {xi) and Pb(xi) for those 
variables. We use the relative fractions in the four lep- 
ton quality categories in the Monte Carlo samples for 
P s (xi) and Pb(xi) for the primary and secondary lepton 
variables. 
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We make a minimum La requirement for each of the 
three signal categories (electron, muon, and pion) which 
has been optimized to give the lowest signal branching 
fraction limit under the assumption of no signal in the 
data. The signal region in m T is defined to be [1.65, 
1.90] GeV/c 2 , which contains 90% of the signal. The 
signal selection efficiency (q), including the Lr require- 
ments, in the m T signal region is 3.17%, 2.04%, and 
2.13% for the electron, muon, and pion channels respec- 
tively. The denominator of these is the same for all 
three and includes all tau decays. We have used a simple 
3-body phase space decay model to generate our signal 
Monte Carlo sample. The systematic uncertainty on e, is 
determined by varying the signal and background PDFs 
for each Because the signal branching fraction is de- 
termined from the ratio of the signal and D^iv yields in 
the data, many systematic uncertainties associated with 
tracking, particle identification, and the B tag reconstruc- 
tion cancel. The amount of background, bi, in the m T sig- 
nal region is estimated from the number of events outside 
the m T signal region (the m T sidebands) in the ranges [0, 
1.65] and [1.9, 3.5] GeV/c 2 and the signal-to-sideband 
ratio from the background Monte Carlo sample. 

The signal branching fraction for each channel (Bi) is 
estimated using the relation 



(ni -bi)/(eiS ), 



(3) 



where n, is the observed number of events in the m T 
signal region, bi is the expected background, is the 
signal efficiency, and Sq is a common sensitivity factor 
given by 



So — 



N 



B 



1 



E tag 
Dpiu 
tag , 



(4) 



where Njj^v, Bd^u, and are the fitted yield, total 

branching fraction, and selection efficiency for the Dp<,v 
control sample and etag^ and ^tag" are the -Btag efficien- 
cies for the signal and D[iv samples respectively. The 
last factor (e^g^/e^g") is determined from the Monte 
Carlo samples and close to one (0.922 ± 0.052) since 
the topology of the events in the signal and D[iv sam- 
ples is very similar. We find Nd^v = 867 ± 52 with 
e Dfn , = 0.345±0.008andBu M1/ = (3.29±0.22)x 10" 3 [3, 
which gives S = (7.0 ± 0.7) x 10 5 . 

The m T signal region in the data was kept blind during 
the development of the analysis, to avoid experimenter's 
bias. After all analysis decisions were made, we "opened 
the box" and found 1, 0, and 2 events in the m T signal 
region for the electron, muon, and pion channels respec- 
tively. These totals are consistent with our expectations 
from background only, which are given in Tabic [T] Dis- 
tributions of m T for the data and for the signal Monte 
Carlo sample are shown in Figure [2] The numbers of 
background events in the m T sidebands are consistent 
with our expectations from the Monte Carlo sample. 




m t (GeV/O 



FIG. 2: Distributions of m T after all selection criteria have 
been applied for the data (points with error bars), back- 
ground Monte Carlo sample (main histogram), and signal 
Monte Carlo sample (inset histogram). The dotted vertical 
lines show the m T signal region [1.65, 1.90] GeV/c 2 . 



The central value for the signal branching fraction is 
B = (0.8 ^2 3) x 10~ 5 , where the uncertainties include 
both statistical and systematic sources. The three chan- 
nels were combined by maximizing the likelihood, which 
is defined as the product of three Poisson probabilities in 
rii, where the mean is given by Equation [3J Uncertain- 
ties on the bi, e,-, and So parameters, which determine the 
Poisson mean for each channel, were included by convolv- 
ing the Poisson PDFs with Gaussians in bi, ej, and So. 
The uncertainties on B correspond to the points where 
the log likelihood drops by 0.5 with respect to the maxi- 
mum log likelihood value. We have verified, with a Monte 
Carlo study, that this maximum likelihood technique is 
unbiased and that the uncertainties are a reasonable es- 
timate of one standard deviation. We find a 90% confi- 
dence level upper limit on the signal branching fraction 
of B < 7.7 x 10~ 5 using the prescription of Feldman and 
Cousins [3] for defining the confidence belt. The central 
value and upper limit on the signal branching fraction 
are both limited by statistical uncertainties. 

In conclusion, we present the first search for the for- 
bidden decay B + — > K + Tfi using 383 million BB pairs 
collected by the BABAR experiment. The observed events 
are consistent with the background-only hypothesis. We 
set an upper limit of B < 7.7 x 10~ 5 on the signal branch- 
ing fraction at 90% confidence level. This result can be 
used to improve the model-independent bound on the 
energy scale of new physics in flavor-changing operators 
reported in [|| from > 2.6 TeV to > 13 TeV. 
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(France), BMBF and DFG (Germany), INFN (Italy), 
FOM (The Netherlands), NFR (Norway), MIST (Rus- 
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TABLE I: The number of events in the m T sidebands, N B b, for the Monte Carlo sample and the data; the ratio of background 
events, inside / outside the m T signal region, BG ratio; the expected number of background events, bi, and number of observed 
data events, m, in the m T signal region; and the signal selection efficiency ti for each of the three channels. The first (second) 
uncertainty on ti is statistical (systematic). 

Channel iV sb (MC) N sh (data) BG ratio bi m u (%) 

electron 5.2 ± 1.3 5 0.10 ±0.05 0.5 ± 0.3 T 3.28 ± 0.13 ± 0.22 

muon 0.7 ±0.5 2 0.30 ±0.15 0.6 ± 0.3 2.09 ± 0.10 ± 0.19 

pion 6.9 ±1.6 14 0.13 ±0.04 1.8 ± 0.6 2 2.18 ± 0.11 ± 0.24 
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